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Abstract 

We are concerned with the stability of formation of large bipolaron in a 3-dimensional (3D) 
crystal. This problem is considered in the framework of nonrelativistic quantum field theory. 
Thus, the Hamiltonian formalism, as Frohlich introduced, is employed to describe the bipolaron. 
We approach the problem by characterizing some sufficient or necessary conditions for the bipolaron 
being stable. This paper gives a full detail of the author's talks at ESI, RIMS, and St. Petersburg 
State Univ. in 2005. 
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I. INTRODUCTION 



The question whether bipolarons can be formed has attracted ever a great deal of studies 
since Shafroth- in particular showed bipolaronic superconductivity takes place when the 
temperature is below that of Bose-Einstein condensation. Over the past few decades there 
have been several renewals^ 1 ^ of interest in this problem from the point of view of high- 
temperature superconductor. In this paper, we will focus our mind on the stability of 
formation of large bipolaron. 

We will treat two electrons coupled with longitudinal optical (LO) phonons in a 3- 
dimensional crystal. Generally, the electron-phonon interaction dresses an electron in a 
phonon cloud. This dressed electron is the so-called polaron. We assume the electron- 
phonon interaction is described by the second quantization, as Frohlich^ introduced. If the 
Coulomb repulsion between the two electrons puts them so far away from each other that 
each electron dresses itself in an individual phonon cloud, there is no exchange of phonons 
between the two. Namely, there is the only Coulomb repulsion between them and thus only 
two separated single polarons are formed in the crystal. On the other hand, if the distance 
between the two electrons is so short that the phonon-exchange occurs, there is a possibility 
that attraction appears^ between them and therefore that they are bound to each other. 
The bound two polarons is called a bipolaron. 

In the light of superconductivity, many studies have been performed for the large 
bipolaron.-^^ -^ 1 ^ As for small polaron,— liiiiLiSiil Alexandrov and Kornilovich^ clarified 
the physical properties of the small Frohlich polaron. Alexandrov and Ranninger- showed 
the possibility that small bipolarons might be superconducting. Moreover, Alexandrov and 
Mott, 19 and Alexandrov 20 pointed out the mobility of the small bipolaron as well as that of 
the large bipolaron. Thus, we are also interested in how the size of large bipolaron grows. 



Using the classical picture described in Ref. 



21 



we can explain the occurrence of attraction 



between the two electrons in the following. Since an electron has negative charge, a first 
electron leaves behind a deformation trail in the crystal lattice, which affects the positions 
of the ion cores. This trail is associated with an increased density of positive charge owing 
to the ion cores. Therefore, it has an attractive effect on a second electron. Namely, the 
lattice deformation causes the attraction between the first and second electrons. Thus, the 
problem is whether the attraction is enough to bind the two electrons. 
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Pekar 22 showed that there is no formation of the large bipolaron for any optic dielectric 
constant eoo and the static dielectric constant e in 3D systems. Also Takada 2 ^ reached 
the same result in 3D systems from the weak-coupling limit by studying the Bethe-Salpeter 



equation. Obeying the classical theory of phonon described in Refj24j, we briefly see that 
two large polarons cannot be bound to each other from the point of view of binding energy 
2Eg™(r) — E^(r), where E^{r) is the total energy of a large (single) polaron in a sphere 
of radius r > 0, and E£^(r) the total energy of a large bipolaron in the sphere. Because 
the binding energy is always negative in classical theory. More precisely, with the help of 
expressions of E£™(r) and E£™(r) given in RefJ24l. we can easily derive the equation: 

- (2E^(r) - Eg{r)) = V2a f-j- - l) ± (1.1) 

for a (the coupling constant of electrons and the phonon field) and 77 := e^/eo (the ionicity 
of the crystal) under the natural units. So, we can say that the binding energy is negative 
in the classical theory because < 77 < 1. A similar form to RHS of Eq. fll.ljl will play an 
important role in Secs JT7| N\ and|Vn](see Eqs.flSD, 0> (E2D> ff57T3|) . and (157141) 1 Thus, 
the attraction derived in classical theory is too weak to bind the two polarons. 

On the other hand, Vinetskii and Giterman 25 found that the bipolaron is formed only for 
very small ionicities 77. Emin 5 also pointed out a possibility of formation of large bipolaron 
for sufficiently large coupling constants a in addition to small ionicities 77. Our argument 
using Frohlich's interaction in this paper is based on many studies in quantum theory by 
prior literature^^ 1 ^ to derive the effective attraction. Bassani et al.— showed the large 
bipolaron formation for large a and small 77 by estimating the binding energy. Verbist et al— 
showed the stable region of 77 for the formation of the large bipolaron in the strong-coupling 
limit. It is worthy of note that the Feynman path- integral approach^ 1 ^ is also useful for 
the bipolaron problem. 

In this paper, to derive an effective attraction in the same framework as in studies,-^ 1 ^ 1 ^ 
we will adopt the notion of Feynman's classical virtual phonon^I instead of the deformation 
trail in the crystal lattice. In addition, we will take the image of Peeters and Devreese's 
classical bipolaron^ into our argument. In their theory of the classical bipolaron, the two 
electrons rotate in a circle around a common fixed center. We now denote the radius of the 
circle by r PD . Then, the classical bipolaron has to be localized in the closed ball B(r PD ) of 
radius r PD centered at the fixed point. The radius r PD should turn out long when the Coulomb 
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repulsion beats the attraction caused by the phonon field. Conversely, it should turn out 
short when the attraction wins over the Coulomb repulsion. So, if the attraction between 
the two electrons is so weak that the Coulomb repulsion makes the interior of the ball B(r PD ) 
completely contains the whole crystal, then we cannot hope that the classical bipolaron is 
formed in the crystal. One of our ideas is to take this image into nonrelativistic quantum 
field theory by reorganizing the way we have ever done.— Thus, we have to introduce the 
notion of the size of the region in which the two electrons constructing bipolaron live. For 
the ground state energy E BP of bipolaron and the ground state energy E SP of single polaron, 
we introduce an energy E(r) coming from the phonon field for every r > as: 

2E SP -E BP = E(r)-j. (1.2) 

That is, we express the binding energy with the energy E(r) and the Coulomb repulsive 
potential U/r. Here, r stands for the distance between the two electrons now. Taking 
account of Emin's work— and Salje's,— E(r) is probably non-negative. Thus, main purpose 
in this paper is to estimate E(r) so that the stability of formation of large bipolaron can be 
characterized in terms of the estimated E(r) and also t]. 

In SecJTIl we will introduce the Hamiltonian of the bipolaron which has the Frohlich 
interaction. In Sec JIIII we will show a device to introduce a parameter 9 > which controls 
the coupling strength between the regimes of the weak- and strong-coupling theories. To find 
this parameter, we will develop Lieb and Thomas' method.— In Sec JIVI we will argue spatial 
localization^^ of the relative motion of bipolaron in the weak-coupling regime by estimating 
E(r). Especially, we will estimate the distance between the two electrons in the bipolaron 
from below by adopting Lieb's idea^ into nonrelativistic quantum field theory. Thus, we 
will understand how the ionicity raises the size of the bipolaron in the weak-coupling regime. 
In SecJVl combining the notion of Feynman's virtual phonon^ and the image of Peeters and 
Devreese's classical bipolaron, 28 - we will derive two effective Hamiltonians in quantum theory 
from the original one. The effective Hamiltonians describe those in the strong-coupling 
regime. With the help of these effective Hamiltonians we will find a sufficient condition for 
the bipolaron formation in terms of E(r) or rj and give a lower and an upper bounds to 
the ground state energy E BP . In Sec JVIl we will argue the spatial localization in the strong- 
coupling regime. In Sec iVIII we will consider a sufficient condition for the positive binding 
energy. 
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II. HAMILTON! ANS 



The total energy of the bipolaron consisting of two electrons coupled with the LO phonons 
is described by the Hamiltonian H BP : 

H BP = H c \- e \ + ifph + -ffel-ph, (2.1) 

where H e \- e \ is the energy of two electrons with the Coulomb repulsion between them, Hp^ 
the free energy of the phonon field, and .ffei-ph Frohlich's interaction^ derived through the 
second quantization of the Coulomb long-range interaction: 

2m J \x\ — X2\ 
H ph = ^ ?^ka{a k , (2.3) 

k 

H el - ph =J2H {Vh^a* + V k *e- ik ^al\. (2.4) 

3=1,2 k ^ ' 

In Eq. (12. 2 p , the position and momentum operators of the jth electron (j = 1, 2) of mass m are 
denoted by Xj and pj, respectively, so pj = —iftV Xj . The strength of the Coulomb repulsion 
is designated by the symbol U, so U = e 2 /e oc for the electric charge e and the optic dielectric 
constant e^. In Eq. fl2.3p . a k and a\ are the annihilation and creation operators, respectively, 
of the LO phonon with the momentum hk. Then a k and a\ satisfy the canonical commutation 
relation, [a^, a\] = 5ke, because phonons are bosons. The LO phonons can be assumed to 
be dispersionless, Uk = uj lo . In Eq. fl2.4p . Vk is given by Vk := —ihuj LO [ATcar fp /k 2 vY^ 2 for 

1 /2 

the crystal volume V and the free polaron radius r fp = (S/2mw LO ) . The dimensionless 
electron-phonon coupling constant is given by 

1 e V 1 1 V (on 

a := — — , (2.5) 

hw LO 2 V e oo eo/ r fp 

where eo is the static dielectric constant. 

Concerning the Hamiltonian H BP of bipolaron, we make some remarks. Since the ionicity 

1] of the crystal is defined by rj := eoo/eo, it satisfies < rj < 1. In terms of r], the strength 

of the Coulomb repulsion is rewritten as U = y/2a/(l —f])- Since we consider the two-body 

system of large polarons, the wave vector k in runs over the first Brillouin zone.— i2£i2I 

This fact makes some noticeable differences between phonon and photon (cf. Table 23.4 of 



Ref 



35l ). The primitive cell is usually given by the first Brillouin zone for a crystal. However, 



it is not always to chose the first Brillouin zone for another solid. When a solid has metallic 
properties, we have to take the Fermi surface into account. In this case, we employ a reduced 
zone scheme.— 

We use the natural units h = m = u ho = 1 from now on. Using the conversion 3 ^ of sums 
to integrals, we estimate \ Vk\ 2 at \^2aK/n as: 

E„ N f o 4:iiar r V f , 47iar r \Fla 
W^^fA dk "777^ = 77T^ / d 3 k — f = — — —K, 2.6 
k V* J mK Vk 2 (2tt) 3 J mK Vk 2 7i 

where V* is the volume of the primitive cell in the reciprocal lattice and V* = (27i) 3 /V c for 
V c , the volume of the primitive cell in the direct lattice. We denoted by N the number of 
primitive cells which is contained in the crystal volume V. So, when integral J d 3 k is over the 
first Brillouin zone, K means the radius of a sphere of the first Brillouin zone. For example, 
in the Debye interpolation scheme 3 ^ for a harmonic crystal, K is given by K = ko, which 
is defined by fc D = (3/47r)(27r) 3 iV/V = 67r 2 N/V. This equation says, as is well known, there 
is the relation between k^ and kp (the radius of the Fermi surface), i.e., kn = (2/Z) 1 / 3 /cf, 
where Z is the nominal valence. When the integral J d 3 k can be extended from the first 
Brillouin zone to the whole fc-space, K plays the role of an ultraviolet (UV) cutoff. 

Using the approximation 3 ^ of the Fourier expansion, V A /(47r|x|) ~ Ylk e x /k 2 , we obtain 

X>|V«^. (2.7) 

We often use this approximation (12.71) in this paper. 

We also use the coordinate of the center-of-mass, X\ = {x\ + x 2 )/2, and the coordinate of 
the relative motion, X 2 = x\ — x%, in this paper. Each momentum is given by Pj = —iVxj, 
j = 1,2. So, we have P\ = pi + j>2 = M\X\ and P 2 = (pi — p 2 )/2 = M 2 X 2 , respectively, 
where we set the masses as M\ = 2 and M 2 = 1/2. Then, H BP is unitary-equivalent to 

H BP = H e i- e i + ifph + -ffel-ph, (2.8) 

where 
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H ehph = ^c(X 2 ,A;){^e M a fc + U;e-* Xl 4}, 
k 
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with 



c(X 2 , k) — 2 cos 



2 



+ e 



-ik-X 2 



As for i? B p, we note the following. The bipolaron Hamiltonian H BP and the total mo- 
mentum Iltot = Pi + Vi + J2k k a k a k are commutable, i.e., [Iltot^Bp] = 0, because of the 
translation invariance of H BP . We can show that H BP has no ground sate in the stan- 
dard mathematical representation because of the continuous symmetry of the translation 
invariance.— It is known in general that if a Hamiltonian H has a continuous symmetry 
(that is, [H, II] = for the generator II of a transformation) and moreover H has a ground 
state, there are two possibilities: the transformation e %m might leave the ground state invari- 
ant for every real number t, which is called manifest symmetry, or otherwise, it might map 
a ground state to another ground state e %m ^o, which is called hidden symmetry. In the 
latter case, we can show that there are infinitely degenerate ground states.— Moreover, the 
Nambu-Goldstone theorem suggests that the spontaneous symmetry breaking would occur 
and then the Nambu-Goldstone bosons would appear. We conjecture that they appear as 
acoustic phonons in bipolaron, taking account of the result in Ref 140 . 

To avoid such a kind of situation coming from the continuous symmetry and thus to give 
a possibility that H BP has a ground state, we can consider, for example, the case where 
we restrict the electrons' movement into the crystal and the case where we nail down the 
center-of-mass of the two electrons at a point Q. For both cases, we only have to employ 



instead of i? e i-ph in H BP , where p(x) is a function satisfying < p(x) < 1. Namely, in those 
cases the total Hamiltonian is H BP = H c \- e \ + H p ^ + /Zei-ph(p)- In the former, p(x) is defined 
so that p{x) = outside the crystal and p(x) = 1 inside the crystal. In the latter, p(x) is 
defined by p(x) = pq(x), where pg{x) :— 1 if x — 2Q; Pq(x) := if x ^ 2Q. Then, we can 
separate the center-of-mass motion from H BP and thus obtain its relative motion. Taking 
p(x) = 1 for every x restores the Hamiltonian i? e i-ei + H ph + H e i- ph (p) to the original H BP . 
When we introduce p(x) into H BP , we employ i? e i-ph(p) instead of if e i- P h, where if e i- P h(p) is 
given by putting p(2Xi) in front of in the definition of i? e i-ph- When we fix the center- 
of-mass at Q, our target Hamiltonian is H BP = H e \- e \ + + H^-^ipo)- As used in Refjol, 
there is another method, which is given by fixing the total momentum n tot at a real number 
in the spectrum of n tot . For example, in the case where Q is the origin O, we can obtain the 



H ehph (p) = p(x 1 + x 2 ) {Vk^ak + V k *e 




} 



3=1,2 k 
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Hamiltonian H%=° := e iP ^H BP e- iP ^ = P 2 /2M 2 + U/\X 2 \+H ph +P 2 h /2M 1 +H chph ( Po ), 
where P ph := J2k ka\a k . 

ILL WEAK- AND STRONG-COUPLING REGIMES 

In this section, a parameter 9 > is introduced to control the strength of the coupling 
constant a in the Hamiltonian of single polaron. We will use this parameter 9 for the 
Hamiltonian of bipolaron as well as that of single polaron throughout this paper. 

Let us define E BP , the ground state energy of bipolaron, as that of H BP , i.e., E BP : = 
inf Spec(H BP ). We denote the energy spectrum of a Hamiltonian H by Spec(H). The 
ground state energy of H always means inf Spec(if) throughout our argument. 

We denote the ground state energy of single polaron by E SP , i.e., E SP : = inf Spec(-ffgjr), 
for the single-polaron Hamiltonian: 

= \$ + E a ^ ak + E { V ^ kXia * + V k*e- lkXj 4} , 3 = 1, 2- (3.1) 

k k 

Through many studies^ 1 ^ 1 ^ 1 ^ 1 ^ we know there are two regimes with respect to powers of 
a in E SP : E SP oc —a in the weak-coupling limit and E SP oc —a 2 in the strong-coupling limit. 
When the wave vector k in Ylk runs over an infinite lattice, Lieb and Thomas^ controlled 
the UV cutoff K by a rigorously to give a lower bound to E SP in the strong-coupling theory. 
But, since our K is now the radius of the first Brillouin zone, we add a device to their idea. 
Using a unitary operator, we introduce the preceding parameter 9 > to control the coupling 
strength between the regimes of the weak- and the strong-coupling theories. Namely, we 
represent the two regimes by 9 in the following. We fix a in the weak-coupling theory first 
(9 = 0). As we switch on 9 (9 > 0) and make it large, we have the electron-phonon coupling 
constant ag := a9 in the strong-coupling theory and ae increases. The method proposed 
here will be basically employed to obtain effective Hamiltonians for the bipolaron in SecJVl 

We define a canonical transformation with the generator G{9) := i9^2, k {Vkdk — Vk a \}- 
Then, our unitary-transformed Hamiltonian is Hq P (9) := e lG ^ H$ P e~ lG ( e \ We note that E SP 
is the ground state energy of H^(9) since Hq P (9) is unitary-equivalent to Hg P . We show 
that E SP can be estimated in the following. It is easy to show that for a in the weak- coupling 
regime 

- < E SP < 0. (3.2) 

7T 
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We can show that E SP oc — a for small coupling constants a in the weak-coupling regime 
by obeying a perturbative method as well known or using a nonperturbative one such as in 
£ Meanwhile, it is shown in this section that for so large 9 that a>g := a>9 



Theorem 2 of Ref . 



43 



is in the strong- coupling regime, a lower and an upper bounds to E SP are given as: 



C sp Q>q 



1 



V2 



K < E SP < - c SP a 2 e + 6 



V2 



Og 



K, 



(3.3) 



7T 7T 

with the approximation [2.1), where c SP = 0.108513- •-. We note that the estimate ( 13.31) 
in the strong-coupling regime tends to the estimate ( 13. 2\\ in the weak-coupling one as 9 
approaches 0. 

To show the inequality (13.31) we need some mathematical arguments. We now introduce 
a mathematical parameter R! > into so that liniR/_»oo Hgp (R 1 ) = (in the norm 
resolvent sensed) as follows: 



H^(R') := \p) + °U + TiM E fa**'"" + V k^' ikXi 4} 



J = 1,2, (3.4) 



where 71 (x) is a smooth real- valued function which satisfies < 71 (x) < 1 and 71 (x) = 
7i(— x) for every x. The parameter i?' > is introduced so that 71 (x) = 1 for \x\ < R'/2; 

for \x\ > R' . We denote by E SP (R') the ground state energy of (R 1 ), i.e., 



7i W 



E SP (R') := inf Spec(H^ (R')) . Then, it automatically follows that liniR/^oo E SP (R') = E SP . 

In the same way as for H$(6) we define H$(6, R') by H$(6, R') := e iG ^ (R')e~ iG ^ . 
Then, through the approximation (12 .TP we approximate H^}(9,R') as: 



p^-29 ll (x J )J2\V k \ 2 e 



+H ph + £ [v k ("A-r ,)<"'■'■ -9)a k + V k * (',(.0), - 9) 4} 

jfc 

We estimate E SP (R') from below in the first step and from above in the next step. After 
obtaining both estimates, taking R' — > 00 yields the estimate (13.31) . 

We arbitrarily fix a normalized phonon state \l/ and a normalized wave function if)(xj) of 
the electron satisfying ]xm\ Xj \-+ 00 'iJ)(xj) = 0. The brackets (( )) stands for an averaging over 
the wave function if)(xj). We define a unitary operator Uj )Q _e so that U* a9 pjUj ta e = a9pj and 
Uj t a6 x 3Uj,a6 = Xj/a9. Then, the term ((U* a0 Hs P (9 , R')Uj, a g}} in the equation, 



(t^V)* HH\9,R') (fW)^ 



((Ul ae Hs P {9, R')Uj t ae)) 



can be approximated as: 

({U* ae HH\e,R')U j!a e)) 
« ^fiipj)) -29j2\V k \ 2 (Mx j /a9)e ik ^ e )) + # 2 ^|^| 2 

fc k 

+H ph + J2{ V d-Yi{xj/ae)e ik ' Xj/ae ~ *K + ^(frifo/a^"*^ - 0)>4}. 
Thus, it is estimated from below as: 

iui ad Hg\e^)u ha8 )) > ^f{(p))) -^iy fe | 2 |((Ti(^/^)e*^})r 

+ (^-l)^|^| 2 |(( 7l (x,/^)e^^))| 2 . 
Since |((7i(x i /a^)e ifc '^/ Qe ))| < 1, we arrive at the inequality: 



> 



({U j>ae iP) *| Hg\e, R') I {U j>a6 iP) *) 
^- J d 3 x\V x ifj(x)\ 2 



|2 ik-(x-y)/a6 



-oJ2\ Vk \ 2 J J d 3 xd 3 yil (x/a9) 11 ( y /a9)\i;(x)\ 2 \i;(y)\ 2 e i 

A; 

-i^-iiEi^i 2 - ( 3 - 5 ) 

Approximation (12 .7p shows 

^|^.| 2 / / rf 3 xA7i(^)7i(y/^)|^W|>(y)| 2 e ifc - ( ^ )/ae 



< -= / / tfztf'y 1 V 1 , 1 , (3.6) 

where inequalities, < ^(x / a9) , ~fi(y / a9) < 1, were used in the above. Thus, the inequali- 
ties, (13. 5p and (13.61) . lead us to the estimate from below: 

((CW) * | H { s i\9, R') | (U j>a6 iP) ) > a 2 9 2 SM) -\9-l\J2 W, 

k 

where S P (ip) is the Pekar functional,— d£dLi& i.e., 

£ P (v>) : = ~y c/ 3 x|v^(x)| 2 - y y 



|^(x)| 2 |^(?/)p 



|x - y| 
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Lieb^ proved that there is a unique and smooth minimizing ip{x) in c SP := — inf^,//^\\ = i S P (ip) 
up to translations. Then, using the estimate (12.61) and the estimate of c SP shown by Miyake 45 
and by Gerlach and Lowen,— we obtain 



Esp(R') > -c SP a 2 e 



1 

1 - T 



X ' /)ae K, c SP = 0.108513- •■ . (3.7) 



7T 



This inequality works for the upper bound in the inequality (13.31) . 

Let us proceed to the proof of the upper bound next. For the estimate from above, we 
prepare 

n<${e, R') ■= \v) - ie lx { Xj ) \v k \ 2 e^ « \v) - ^ + 

k 

V2a9 



\Xj | 



error ^(x^tf), 



eTTOT Rl (xj,9) : = (1 - Jiixj)) 



\Xj I 



by using the approximation (12.71) . 

For r > we define a function Cq(t) by 

C (r):=(l + r)e- r . (3.8) 

Since C (r') < C (r) < C (0) = 1 for < r < r' and C (oo) = lim^^ C (r) = 0, we take a 
unique point r* so that C (r*) = 1 — v / c^'. We set Ci(/i, i?') and as Ci(/z, i?') := C (fiR') 
and := r*///, respectively. Then, we have a simple inequality: 

d^i?')^!"^ (3-9) 



provided R^ < R' . Using the wave function ip^Xj) defined by ip^Xj) := a/ /i 3 /vr e At ' x -''' for 
/i > as a trial function, we can estimate (ipJ H^(9, R') from above: 

(W| H^{9,R') \%) < ^ - V2a ef x+ v^o^C^tf), (3.10) 

where we used < 1 — 7i(x,) < 1. Evaluating RHS of the inequality E SP (R') < 
(Vv^phl Hsp{9, R') |^fiph) for the phonon vacuum f2 p h, it follows from the inequalities 
(13791) and (15301 that 



E SP {R') < \ (fi-V2^ae) 2 -c SP a 2 9 + 9 2 J2\Vk\ 
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for every p > and R' > R^. Setting the parameter p as p = /i* := y/2c SP ae and using 
estimate (12. 6p . we arrive at the following upper bound: 

E SP (R') < inf (^fiphl #$(0, iO |V^ph) < -c SP « e 2 + (3.11) 

/i>0 7T 

Both bounds (13. 7p and (13. lip lead us to the conclusion that if R' > R p>t = r*/p* = 



r*/ V2c SP «e and > then 

1 



^ ae K< E SP (R>) < -c SP a 2 d + 9^K. (3.12) 



7T 7T 



Taking the limit R' — ► oo in the above, we reach our desired estimate (13.31) . 

IV. SPATIAL LOCALIZATION IN WEAK-COUPLING REGIME 

In this section, we consider the case where a is in the weak-coupling regime. As noted in 
SecJTTl the Hamiltonian formalism in the mathematical representation does not show us the 
existence of any ground state in the state space because of the translation invariance of H BP 
with p(x) = 1. Thus, we employ H e \- ph (p) in H BP so that p breaks the translation invariance. 
Then, in the weak-coupling regime, we seek a necessary condition for H BP having a bound 
state and investigate spatial localization of bipolaron. We note that we can obtain the same 
results for H™ p l and H^ as described below. Although we do not prove this remark in this 
paper, it is proved in the same way as shown for H BP . 



In Refj29l we showed that if the quantum particle is not spatially localized, then the 
particle dressed in the cloud of bosons cannot exist. In that problem, divergent soft bosons 
for the infrared catastrophe cause the spatial nonlocalization. From this point of view, we 
study how the Coulomb repulsion between the two electrons arrests the spatial localization 
in bipolaron. 

Let us proceed with the definition of spatial localization. We say that the rela- 
tive motion of the bipolaron in a bound state \l/ n is spatially localized in the closed ball 
B{r) := {X2 I IX2I < r} if (\l/ n |\]/ n ) 1 (\l/ n ||xi — a^lj^n) < t. In this paper we will use the 
phrase that a bound state \l/ n of the bipolaron with the size less than r > does not exist 
when the relative motion of the bipolaron in \l/ n is not spatially localized in B(r). For a 
ground state \l/ we define the distance d BP (^ ) between the two electrons in bipolaron by 

d BP (^ ) ■= <^ |^o) _1 <^o|ki -ar 2 ||*o>- (4-1) 
12 



To estimate E(r) in Eg. (11.21) under the weak-coupling regime, we define a positive con- 
stant E w (a) by 

Bw(a): =4£|n| 2 =^K=^ (4.2) 

for every a > 0, where A is a wave length defined by A := 2ir/K. It is easy to check 
E BP < by evaluating RHS of the inequality E BP < inf^ ^/>f2 ph | H BP \ipfl ph ) for normalized 
wave functions ip of the two electrons and the phonon vacuum f2 p h- The simple operator 
inequality H pil + H e i- ph > — E w (a) shows that H BP > — E w (a), making a lower bound to 
the ground state energy: 

E BP >-E w (a). (4.3) 

We can show E BP oc — a for small a in the weak-coupling regime by a perturbative or 
nonperturbative means in the same way as in the case of E SP . 

We prove in this section that if the bipolaron has a bound state with eigenenergy E so 
that its relative motion is spatially localized in B(r), then there is the energy inequality: 

E w (a) + E>-= V " . (4.4) 
r (1 — r])r 



Our result says that 



d BP (*o) > C := F ( f, F > ^7T^ (4-5) 
E w (a) + E BP 8(1 - rj) 



for any ground state ^>q. Thus, we realize that the distance d B p(^o) between the two electrons 
in the bipolaron grows more and more as 1] turns out to be closer to 1 and also that the 
coupling constant a in the weak-coupling regime cannot stem its growth. That is not the 
case in the strong-coupling regime (see Sec lVIl) . This is a noticeable difference between the 
cases in the weak-coupling regime and in the strong-coupling one. 

Moreover, this growth tells us that if a bipolaron has a ground state for rj sufficiently 
close to 1, then B(d BP (^fo)) completely contains the whole crystal. On the other hand, since 
the bipolaron in the ground state should be formed in the crystal, the two electrons 
in the bipolaron must exist in the crystal. Thus, we have met a contradiction. Therefore, 
reductio ad absurdum makes us conclude from these two facts contradicting each other that 
any bipolaron in a ground state cannot be formed in the crystal if r\ is close to 1 in the 
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weak- coupling regime. We will give another lower bound to rf B p(^o) i n the inequality (I4.10P 
later. 

Our above result also states, in particular, that the relative motion of the bipolaron in 
a ground state is not spatially localized in B{r), provided that the ground state exists under 
the condition: 

E w {a) < U - (i.e., < 1 - A < ^ . (4.6) 

In terms of the binding energy with Eq. fll.2l) . the condition (14.61) means 

- (2E SP - -£/ B p) « - E w (a) -E BP + - 

r 

> -E w (a) + - = V2a(-^-^-)->0 (4.7) 
r \1 — rj A J r 

by estimating E(r) at E w (a) + E BP . Thus, as for the relation between the spatial localization 

of a ground state and the binding energy, we can interpret our result as that a ground state 

of the bipolaron with the size of d BP {^o) ^ ess than r satisfying the condition ^4.6)) cannot 

exist because of the negative binding energy fl^.Tp . 

We define the number of phonons in a ground state \l/o by N p ^(^fo) := 

(\l/o|^o) 1 (^ f o| Sfc a l Q fc|^o)- Then, it is shown in this section that iV ph (\I/ ) is estimated as: 

, 7^ . + (2E SP - E BP ) < N ph (^ ) < 2E w (a) + E BP + 2^ E w {a){E^{a) + E BP ). (4.8) 

"bpI^oJ 

The first inequality yields the upper bound to the binding energy: 

2E SP - E BP < iV ph (^ ) - , Wl : < iV ph ($ ). (4.9) 

"bp^oJ 

Hence it follows that another lower bound to d BP (^ ) is given as: 

wrk^i £ M (410) 

Let us begin by proving the first result. We adopt Lieb's idea^ into our argument to 
prove it. Here we assume that H BP has a bound state \l/ n with its eigenenergy E, i.e., 
H BP ^ n = E^f n , and that it is spatially localized in B(r). Here we can suppose \l/ n is 
normalized, i.e., (^f n r& n ) = 1) without any loss of generality. Because we only have to take 
i^! n ^n) as the normalized bound state in case that \l/ n is not normalized. Then, 

since P^/2M 1 > 0, evaluating RHS of = (\X 2 \^ n \H BP - E\^ n ) from below yields 

> A-(|X 2 |* n |P 2 2 * n ) + U + {\X 2 \ l ' 2 ^ n \H p)x + H^WX^'^n) - E(* n \ \X 2 \ |* n >. 



2 
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This leads us to the inequality: 



> ^-<|X 2 |vI/ n |P 2 2 vl/ n ) + U-(E + P w (a)) r, (4.11) 



2 



where we used the spatial localization of and the fact that E + E w (a) > E BP + E w (a) > 
by the inequality (14. 3p . It follows from the same argument that 



X 2 \V n ) + U -(E + E w {a))r. (4.12) 



> -^(P^n 
~ 2M 2 \ 2 

Since &x 2 \X 2 \- 1 = 4tt5(X 2 ), we have (Pi\X 2 \ + |X 2 |P 2 2 ) /2 = \X 2 \P 2 \X 2 \- 1 P 2 \X 2 \ - 
27r|X 2 |5(X 2 )|X 2 |. Thus, we reach the equation and inequality: 



P 9 I Xo I + I Xo I Po 



P2\X 2 \V n ) > 0. (4.13) 



2M 2 

Combining the inequalities (I4.11I) -( T4.13I) . we arrive at the first desired inequality (14.41) . 

To show the inequality (14.51) . we use reductio ad absurdum. As an assumption for it we 
suppose d BP < d l °p. Then, there is a positive number r' so that d BP < r' < d l BP . Since 
d BP < r', the above fact already proved says E w (a) + E BP > U/r', and so r' > d x BP . This is 
not consistent with the fact that r' < d l BP . Thus, reductio ad absurdum supplies the second 
desired inequality. 

Let us proceed with the proof of our next statement. Since E BP < 0, the simple term 
E w (a) can substitute for E w (a) + E BP in the inequality (I4.4p for the ground state ^/ - 
Thus, the contraposition of the first result leads us to the third desired result on the spatial 
nonlocalization. 

Let us estimate the number of phonons in a ground state now. We regard as 
being normalized. We define the Hamiltonian 52j Hs P of the two separate single polarons 
for H BP by using the coordinates X\ and X 2 , i.e., 52j Hsp := J2j=i 2 -P//2M 2 + 2P p h + -£/ei-ph- 
In addition, we note that the Hamiltonian of the two separate single polarons for H BP is 
Efifs? = Pf/2M 2 + 2P ph + Pei-ph (p Q ), and that for P B n =° is 5^ n=0 #sp = Pl/2M 2 + 
2PT ph + P 2 h /2M! + PTei-ph(po). 

Note 1 = |X 2 | _1 / 2 |X 2 | 1 / 2 now. Then, Schwarz's inequality implies that 1 = 
(iXal- 1 / 2 ^)!^! 1 / 2 ^) < (^ollXal- 1 !^) 17 ^^!!^!!^) 172 , making the inequality 
l/d BP (^>o) < (^oj |X 2 | _1 |\l/o)- Hence it follows from this inequality, together with sand- 
wiching the equation 52,- Hsp ~ 2-^sp = -£^ B p — 2P SP + if ph — U /\X 2 \ between (^o\ and j^o); 
that < (P BP — 2P SP ) + iVph^o) — U / d BP {^! 0) . Therefore, the lower bound in the inequality 
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(14. 8 p is obtained from this inequality. As for the upper bound, we use the simple operator 
inequality: 

CH BP - H ph > (C - 1) j# ph + G _ ^ Fei-phj > - -^—jE w (a), C > 1, 

to come up with H p ^ < CH BP + C 2 E w (a)/(C — 1). Putting this inequality between (\I>o| 
and |^o) leads us to the inequality N ph (^ Q ) < inf c >i [CE BP + C 2 E w (a)/(C - 1)]. RHS of 
this inequality attains the lower bound in the inequality (14.81) . 



V. BIPOLARON FORMATION 



In this section we deal with the original H BP (i.e., in the case p(x) = 1) in the strong- 
coupling regime and derive two effective Hamiltonians from H BP by modifying Bogolubovs' 
method.— Our modified method is similar to Adamowski's^ and ours.— More precisely, 
we seek a canonical transformation Uq with the parameter 9 > so that the transformed 
Hamiltonian has the form of H BP (9) := UgH BP Ug = H ef[ (9) + H ph + H c \- ph (9) + where 
H c fi(9) is an effective Hamiltonian in quantum mechanics, and a divergent energy as 
9 — » oo. The canonical transformation Uq requires that the effective Hamiltonian H e g(9) 
should gain an attractive potential V(9) from the phonon field as H c r(9) = H c \- C \ + V(9). 
Thanks to this extra attractive potential V(9), we expect that there is a critical point 9 C 
so that the Hamiltonian H e ^{9) itself or the Hamiltonian H T e ^{9) for the relative motion of 
H c q(9) has a ground state if 9 > 9 C . On the other hand, it has no ground state if 9 < 9 C . 
Thus, we can expect that the bipolaron is stable at/near the point 9 C . Actually, according 
to the recent result,— if we apply the above method to -f^T , then we might be able to show 
that H^ has a ground state for sufficiently large 9 > 0. 

To find the canonical transformation Ug, we adopt the notion of Feynman's virtual 
phonon^ and the image of Peeters and Devreese's classical bipolaron,— being led to the 
classical images of the two kinds of states of bipolaron: balanced state and unbalanced state. 
The balanced state designates the image that the virtual phonon sits at the center of two 
electrons. On the other hand, the unbalanced state represents the image that the virtual 
phonon is not on the line segment between the two electrons. Then, two effective Hamilto- 
nians are derived in quantum mechanics through these classical images. 

It is shown in this section that the approximation \2. 7| ) yields an effective Hamiltonian 
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describing balanced state in quantum mechanics: 

H cS (9) = H ehel + V{9) = \ V \ + \y\ + | U{9) | (5.1) 

2 2 \xi — x 2 \ 

with 

(5.2) 



U(6) = U- V2a9 = V2a (r~ ~ 9 J ■ 



It is clear that 8 C = 1/(1 — rj). Let us fix c* arbitrarily so that < c* < 1. We set c BP as: 



2 / 1 x J 



Cbp-^-^-^J . (5.3) 

Then, using the effective Hamiltonian H e ^{9) in Eq. fl5.lj) . we show that an upper bound to 
the ground state energy E BP given as: 

2 , V2a ( A ^ , K 3 \ 



provided that 



E,p < - c BP a„ + :L -^ L UK + —\, (5.4) 



1 - ± > „. (5.5) 

We note that the condition (15.51) prohibits us from taking the limit 9 — > in the inequality 
f!5.4p because 

e>^-> 1 



C*(l — 77) 1-7?' 

Meanwhile, it is also shown in this section that a lower bound to E BP as: 

E BP > ( inf Sg(<p)) a 2 e -\l-h ^^-K, (5.6) 

V <p J 1 9 1 7T 



where 



S e {(p) := - / / d 3 Xl d 3 x 2 




Ivy /- \ | 2 , ivy / \|2 . 2v/2 xi,x 2 

^(l — 77) \x 1 -x 2 \ 



~ I j d 3 x l0 l 3 x 2 I I d 3 Vl d 3 y 2 



|^(^1j»2)|' 













(5.7) 



,2 

an energy functional describing unbalanced state and inf^, £o(<p) < 0. This functional is 
the almost same as Vinetskii and Giterman's.— i2^^2i^i Namely, for wave functions (p(xi,x 2 ) 
satisfying (p(xi, x 2 ) = ^>(x 2 , Xi), our energy functional becomes theirs. Here, we note that the 
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estimate (15. 6p in the strong-coupling regime recovers the estimate f!4.3[) in the weak-coupling 
one as 9 = 0. 

At the end of this section we show that the approximation \2. 1\ ) yields another effective 
Hamiltonian describing unbalanced state in quantum mechanics: 



„-_1 9 A— -l 9 I J < \ 



V2 



9(1 - rj)\x\ - x 2 \ 



(5.8) 



U=l,2 3=1,2 

Here we note that the approximation (12.71) breaks the translation invariance so that the 
effective Hamiltonian H e $(9) in Eq. (j5.8j) describes that the virtual phonon is nailed down 
at the origin. Let us set E s (a) as: 

E s (a) := - i - 1 ) a9 (5.9) 

now. Then, the effective Hamiltonian H e g(9) in Eg. / 1 5. 8\) has a ground state provided that 
there is an r > so that 

-<E s (a). (5.10) 
r 

Thus, we are led to the conclusion that 9 C < 1/(2 — v / 2)(l — tj). The condition ( 15. lOf) puts 
restrictions on 9,rj and r. Namely, the sufficient condition for the inequality $5.10}) is: 

1- 1 -^>V, (5.H) 
V29 

and 



Re, v - \]Re, v 2 - K r < Ro tn + ^R diV 2 - 1, (5.12) 

where 

Therefore, for sufficiently large 9 (i.e., 9 ~ oo), 0.585 < r/r fp < 3.415. 

Since the inequalities (13. 3ft . (15.41) and (15.61) say that E(r) > is dominated from below 
by \^2ae/r and/or E s (a) for sufficiently large 9 > 0, in the case of the balanced state the 
condition 9 > 1/(1 — 77) makes the estimated binding energy positive: 

- (2E SP - E BP ) < + E = ^2 a (J— _ 0] I < 0. (5.13) 

r r \1 — 77 / r 
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In the case of the unbalanced state the condition (15.101) also makes the estimated binding 
energy positive: 

- (2E SP - E BP ) < -E s (a) + - = a e 

r 

Thus, the appearance of 9 in Eqs.f l5.13j) and (15.141) is very different than it does not appear 
in Eqs.flU]) and fTCTj) . 



rj) ) r 



< 0. (5.14) 



A. Strategy for Effective Hamiltonians 

In Ref J33| we derived the Coulomb attractive potential from the electron-boson interaction 
for the model introduced by Gross.— Similarly, we use such a canonical transformation with 
trial functions (3j(k,Xj), j = 1,2, which have the parameter 6. By controlling the trial 
functions, we derive the effective Hamiltonian H e s(9). Set (3j^ as j3j t k '■= fij{k,Xj), where 
we assume fij(k,Xj)* = /3j(—k,Xj) for every Xj and j = 1,2, i.e., f3* k = Pj-k- We give 
the generator T by T := i J2j=i,2 Y.k{^Pj,ke ik ' X:i ak - V k * /3* <k e~ ik ' Xj al}. Then, we obtain the 
canonical transformation Ug by Ug := e lT . Precisely writing down, we obtain each term in 
the following. The effective potential V(8) is given by 



V(6) := 2]T \V k \ 2 Ul k ^ k -(#,* + fc,k 

k I 



The electron-phonon interaction H e \- P h(6) is decomposed into the two parts as H e \- P h(9) 
Hehph + ^ei-ph) which are respectively defined by 



-"el-ph • 



E E I Vk f 1 - - \ D elkX3ak 

(=1,2 k I ^ ' 

+V k * (l - /3* fe - \d e 

:= E {pa + A h} + \ E { A ) + 2A] A + A ) 2 } 



3=1,2 3=1,2 



where D(p j)k ) = -A Xj P jtk - 2/7/V,. >',./, + fc 2 ^ and Aj = Y,k Vk(kPj,k - '(V,-, ^,k)je ikx ^a k . 
The divergent energy as # — > oo is 

s * = E l^l 2 | ( x + y) OW 2 + IW a ) + \ (|v,A*r + |v, 2 A fc | 2 ) 

-&(J3i,k + M ?■ 
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The main problem in our strategy is what is the principle to find the best effective 
Hamiltonian. Here we check well-known method for a while by way of trial. If we employ 
(3j(k) = 9(1 + k 2 /2)~ 1 for 9 > (J — 1, 2), the e tT as 9 = 1 is the very unitary operator we 
usee— for the model introduced by Gross. That is, it is the unitary operator in Tomonaga's 
intermediate coupling approximation^ rearranged in Lee, Low, and Pines' study^ 1 ^ of the 
polaron. In the intermediate coupling approximation, the variational principle works so that 
the self-energy turns out lowest. Thus, employing this fij^ = with the help of the 

approximation of the Fourier expansion (12.71) we obtain 

4y/2wa \ fir 0(2 - 0) „,„ m /5F e -^»-»l , „ 2 /7I __^ |ll _ Ja| 



V(6) = " TWf 1 V 2 R^J " 9(2 " 6) V 2 -W^A + ° V 2 ' 
Since £/ = \^2a/(l — t]), the effect Hamiltonian H e s(9) can be derived as: 



1 2 1 2 / 1 „,„ „A V2 



SeffW = rPl + o^2 + 1 0(2 - 



! o 



2 2 \1 — 77 / l^x — X2I 

+#(2 - 9)V2a- - e 2 V2ae-^ xl - X2 K 

Fi - rr 2 | 

Unfortunately, we cannot gain enough attraction from the electron-phonon interaction be- 
cause of the same reason as in Eq. fll.ip . Thus, we have now lost such a principle to derive the 
best effective Hamiltonian as in Lee, Low, and Pine's study. Then, we depend on the picture 
given by the combination of the two notions: Feynman's of virtual phonon, 27 and Peeters 
and Devreese's of the classical bipolaron.— As is in Kadanoff 's explanation,— Feynman con- 
sidered that the electron in the phonon cloud, which forms a single polaron, is coupled to 
another particle with a harmonic oscillator. We call the particle a (classical) virtual phonon. 
We use the notion of this virtual phonon. Although Feynman assumed harmonic oscillator 
for the interaction between the electron and the virtual phonon, we assume the Coulomb 
attractive potential between each electron and the virtual phonon instead as follows: We 
imagine the Peeters and Devreese's classical bipolaron first. We assume a (classical) virtual 
phonon sits at the center (x\ +X2)/2 of the line segment of x\ and x-i. In this case, the two 
electrons can feel an attraction between themselves besides the Coulomb repulsion without 
noticing the existence of the virtual phonon, though the attraction is actually made by the 
virtual phonon. Thus, we reach the image of the balanced state. On the other hand, for 
the unbalanced state, the two electrons have to become aware of the existence of the virtual 
phonon. 
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B. Effective Hamiltonian for Balanced State 



In this subsection, we consider how we can obtain the Coulomb attraction between the 
two electrons in quantum theory for the balanced state. To derive such the Coulomb at- 
traction, we employ f3j^ simply defined by @i t k = —f3 2 ^ = V9- Then, is approximated to 
2 v / 2a(# - \f6)K/-K + V2a6K 3 /67t by the estimate fj2~Bj) . so that lim^*, £ e /# 2 = 0. Then, 
the approximation of the Fourier expansion (12.71) leads us to the more effective V(9): 

\Xi - X 2 \ 

Therefore, we obtain the effective Hamiltonian H e g(9) defined in Eq. fl5.ip representing the 
balanced state. Using coordinates of the center-of-mass and the relative motion, we realize 
that H e s(9) is unitary-equivalent to H e s(9): 

H cS (9) :=\P? + H$(9), 



rel^ ... D 2 , U(9) 



where 

KM := n + , v , ^ 

At the end of this subsection, we show the upper bound (15.41) to E BP . We now introduce 
a mathematical parameter R > 2 into H BP so that lim^oo H BP (R) = H BP (in the norm 
resolvent sensed). To introduce the parameter R > 2, we use a real- valued function 72 (a?) 
defined by by 72(2) := X-r(M); where Xn( r ) '■= 1 for r < R/2 and Xr(t) '■= for r > R with 
linear interpolation. Then, it is easy to show that sup x |V x 72(x)| < 2/R < 1. We define 
the Hamiltonian H BP (R) with the parameter R > 2 by H BP (R) := i? e i-ei + -f^ph + -^ei-ph(-R), 
where H e \- ph (R) is defined by putting 72 (#1 — x 2 ) in front of ^,=1 2 2^2 k ^ n ^ ne original -^ei-ph 
of H BP . 

We can define the unitary-transformed Hamiltonian H BP (9,R) by H BP (9,R) := 
e iT H BP (R)e- iT . Then, H BP (9,R) is approximated as H BP (9,R) « H cS (9,R) + F ph + 
H e \- p h{9,R) + through the approximation ( 12. 7ft . Here the effective Hamiltonian 

H c{i (9, R) is given by H e s(9, R) = H e{i (9) + error H (xi — x 2 , 9) for the effective Hamiltonian 
H e ff(9) in Eq. fj5.lj) and the error term error fl (a;i — x 2 , 9): 

V2a9 



error^zi - x 2 , 9) = (l - ^y 2 {xi - x 2 ) 



Fl _ X 2\ 
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and the electron-phonon interaction if e i_ p h(0, R) is given by //ei-ph(0j R) = H^] ph (R) + H^] ph 
for If^-ph given in SecJVA]and H^ ph {R): 



.7=1,2 fc I V 



" - x 2 ) - - (Pit) ) ' '''''''"> 



\ ] ( 72(^1 - x a ) - fa - \D fe)* ) e-^a[ 



The self-energy T*e(R) in H BP (9,R) is given by 

= l^l 2 | (1 + y) (IW 2 + l/M 2 ) + \ {\VM 2 + |v, 2 A fc | 2 ) 

-72O1 - 2:2)^ ((3i,k + /3 2 ,fe) 

Using coordinates of the center-of-mass and the relative motion, we realize that H e d(9,R) 
is unitary-equivalent to H cS (9,R) := P* /4 + H T J(9, R), where H$(6,R) := # c r ff(#) + 
error H (X 2 , 0). 

For the real number c#, we can take a unique point r BP so that Co(r BP ) = 1 — c*, where 
C (r) was defined in Eq. fl3.8p . We set C 2 {n,R) and i? 

bp as C 2 {n, R) '■ — Co(//i2) and i? BP : — 
r BP //i, respectively. Then C^ii.R) < 1 — c*, provided i? M < i?. Setting r as r = in 



the definition of ip^Xj), we have the wave function M (Xx) as m (Xl) := ^/Jj^Jtt e ^l Xl L 



We denote that of X 2 by ip fl (X 2 ) := y^JjFJn e ^ x ^. Set _R BP as _R BP := r BP //x. Then 
C 2 ([i, i?) < 1 — c*, provided i? BP < -R. It follows from simple estimates that 

for every which satisfies i? BP < i?. It is easily shown that (0j Pf |</> M ) = [J 2 - We define 
a wave function <^(X l5 X 2 ) by ip(Xx,X 2 ) := (j) ^{Xxjij) ^{X 2 ) . Then, for every \i > and 
-R > -Rbp, the ground state energy E BP (R) = inf Spec(H BP (9, R)) of H BP (9,R) is estimated 
from above as: 

£bpOR) < <^ph| H BP (9,R) \ V n ph ) < ^<0 M | Pi + (^| S$(0,.R) M + MR) 



2 ■ 2 



+- — 
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Setting the parameter fi as \i = /x BP := (2<\Z2/5)aO(c* — 1/0(1 — 77)) > because of the 
guarantee (15.51) for it, we arrive at 

E B AR)<-c BP al + ^(4K+!f) 

for every R > max{2, i? BP }. Taking R — > oo in the above, we obtain our desired upper 
bound. 

C. Effective Hamiltonian for Unbalanced State 

In this subsection, we consider how we can obtain the Coulomb attraction between each 
electron and virtual phonon in quantum theory for the unbalanced state. 

Let us proceed with deriving the energy functional £e((f) in Eq. (15.71) first. We arbitrarily 
fix a normalized phonon state \l/ and a normalized wave function ipipc^x^) of the electron 
satisfying \\m\ Xj \^ OQ ip(xi 1 X2) = 0. We denote by ( ) an averaging over the wave function 
(p(xi,X2). Set U a e := U\ t aeU2, a e- Then, in the same way as in the case of single polaron, we 
can estimate the term (U^ 9 H BP U a g) in the equation, 

((U ae tp)y\H BP \(U ae tp)*) = (y 

from below as: 

(K e H BP u ae ) > «pf> + (pi)) + (r-^) -oJ2\v k 

L \X\ — X2\ 

- \9-l\E w {a). 

Using the approximation (12. 7p . in the same way as for E SP , we reach the lower bound (15.61) . 

To obtain an effective Hamiltonian for the unbalanced state, we only have to employ 
G{6) as T. Namely, set f3 j)k as /3 jtk = (3j{k) = e'^O. Then, we have T = G{6). Under the 
approximation (12. 7p . we approximate H BP [9] := U* 9 e lG ^ H BP e~ lG ^U a e to: 

H BP [6] ta H cS (8) + H ph 

+ J2J2{ Vk {e ik ^ /ad -e)a k + V k * ( e - ik ^ /ad - 6) 4} + E e , (5.16) 
i=i,2 k 

where H e s(9) is given by Eq. (I5.8j) and Eq := 6E w (a)/A. We note again that the approxima- 
tion (I2.7P breaks the translation invariance in the original Hamiltonian H BP . 
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(U* e H BP 



jik-xi/aB 



_|_ e ik-X2/ae\ 



Let us fix r > arbitrarily. We denote by ipj(xj) the wave functions a/ (1/irr 3 ) e~\ Xj ^ r 
for j = 1,2. We define the wave functions i/; r (xi,X2) with the parameter r by if) r (xi,X2) := 
ipi(xi)ip2{ x 2) ■ We decompose H eS (9) into the following: 

V2 



H eS (9) 



1 2 

2 Pl 



1 2 



V2 



9(1 - 7]) \\xi - x 2 



1 
1 



1 



F2 



It is well known (see page 89 of RefiSfl) that (ip r \ \ X\ — X2 | ^ — \x 2 I 1 |Vv) < is obtained since 
ipi(xi) is spherically symmetric. Thus, the inequality (ip r \ ag 2 H c g(9) \ip r ) < \/2(l/0(l— ??) — 
2)/r + 1/r 2 < —1 is obtained from the above decomposition and the inequality. Therefore, it 
follows from the HVZ theorem 56 that H e ^(9) has a ground state under the condition (15.101) . 
The restrictions on 8, r], and r are obtained by solving the inequality (I5.10p with respect to 
the variable r. 



VI. SPATIAL LOCALIZATION IN STRONG-COUPLING REGIME 

In this section we consider the approximated H BP [9] given in Eq. (15.161) . For a ground 
state \l/o °f H BP [6], we define the radius u bp (^/q) of the sphere in which the two electron lives 
by 

it BP 

(\[>o) := max 

U#o|*o) (#o|N|#on. (6.1) 

Then, we show in this section that if the bipolaron has a ground state ^o, then there is a 
relation: 

M >2Jl+(l + 5f) ( 1 - 2] . (6.2, 



V^l V 4; a 2 g J \8(1- V ) 
Thus, we are led to the conclusion that, even if 77 approaches 1, 9 > in the strong-coupling 
regime works to stem growth of the u BP (^o)- This is a noticeable difference from the case 
of the weak-coupling regime (14.51) . 

We can show the above result in the same way as we did in Sec JIVl Adopting Lieb's 
idear^ into our argument lead us the inequality, 



V2- 



0>-«kp(*o)- v / 2«2 + -X e - Uj -(4 + 39) 

0(1 - 77) |xi - x 2 \ 
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Here we used the inequalities E BP < andp|,/2 — \/2/\xf\ > — lforj'^j. Since \x\— x 2 \ < 
\ x i \ + 1^1 , we eventually obtain 

> - ^aj + (2 + £ w (a)j m bp (^o) + V^ag " 2 ) ' 

which implies the inequality (|6.2p . 

VII. POSITIVE BINDING ENERGY 

Combining the inequalities (13.31) and (15. 4p . we obtain a sufficient condition for the binding 
energy being positive. Namely, if c* ; 9, and rj satisfy c BP > 2c SP , £/ien the binding energy is 
positive, i.e., 

E BP < 2E SP (7.1) 

for sufficiently large 9 > with the condition A5.5\) . Here, remember c SP = 0.108513 
Then, we note that < c BP < 0.4 and lim c ^i 5 e^oo = 0.4 under the condition 1 15. 5ft . More- 
over, the condition c BP > 2c SP is equivalent to 

5\ 0(1 -V) 



\ ( c * ~ 7777^ — \ ) > 0.217024, (7.2) 



which has the form with corrections in the inequality (12) of Refi5ll. Thus, the condition 
(17. 2p is almost equivalent to 

' > c *~ 77777^75 >V- (7-3) 



0.26350 0.26350 

As Smondyrev and Devreese^- pointed out, we must keep it in mind that the inequality (17. 2p 
is a sufficient condition for the positive binding energy (17. ip . not a necessary condition. If 
we take the effect of not only the leading terms of c SP and c BP but also their remainders into 
the condition (17. 3p . we should control c* so that rj and 9 meet to the results by experiments. 

According to the recent result of study 49 of the Hamiltonians ^sp an d H^ , 

we might be able to choose — inf v ^)) =1 S e ((f) as c BP (i.e., — c BP — inf^ £q{^pY) so that 
2E SP — E BP = — (2c SP — c BP )aj > for sufficiently large > 0. We note that this type of 
equation expressing the binding energy is pointed out by Vinetskii and Giterman.— 2k 
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